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Abstract 

We study a nonlocal aggregation equation with degenerate diffusion, set in a periodic domain. 
This equation represents the generalization to m > 1 of the McKean-Vlasov equation where 
here the "diffusive" portion of the dynamics are governed by Porous medium self-interactions. 
We focus primarily on m £ (1,2] with particular emphasis on m = 2. In general, we establish 
regularity properties and, for small interaction, exponential decay to the uniform stationary 
solution. For m = 2, we obtain essentially sharp results on the rate of decay for the entire 
regime up to the (sharp) transitional value of the interaction parameter. 

1 Introduction 

In this paper we study weak solutions of the equation: 

=A(p™) + 0L'^(2-'")V-(pV(F*p)) inTix[0,oo), (1.1) 

where * stands for convolution, and the space domain is the d-dimension torus with scale L, defined 
as := [— ■j]'' with periodic boundary condition. We assume that V smooth and integrable 
(for precise conditions, see (V1)-(V2) in Section|3]), and that 6* is a positive constant. The primary 
focus of this work concerns the cases m G (1,2] - especially m = 2. In addition, we remark that a 
goal of interest (not always achieved) is to acquire results uniform in L for L 3> 1. We point out 
that, in the absence of the aggregation term (i.e., when V ~ 0) our equation becomes the well-known 
Porous medium equation (PME): 

Pt - A(p'") = 0. 

Note that, formally (and in actuality) the mass of the solution to Ea. (|l.l[) is preserved over time: 

J p{x,0)dx = J p{x,t)dx for all t> 0. 

Without loss of generality, we can thus assume J p{x,0)dx — 1 and results for other normalizations 
can be obtained by scaling. 

In the context of biological aggregation, p represents the population density which locally disperses 
by the diffusion term, while V is the sensing (interaction) kernel that models the long-range attrac- 
tion; Eq. (|l.ip is relevant for models which have been introduced by |BCM| and |TBL| . and further 
studied by |BCM2j and |BF| . The above equation can also be regarded as the evolution equation for 
a strongly interacting fluid: The V represents the long distance component of the interaction while 
short distance interactions - and entropic effects - are accounted for by the degeneracy (m > 1) 
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in the diffusion term. Mathematically, the equation exhibits an interesting competition between 
degenerate diffusion and nonlocal aggregation. 

When V satisfies V{x) = V{—x), Eq. ljl.ip is a gradient flow of the following energy with respect 
to the Wasserstein metric: 

Mp) ■■= I - + * P)dx. (1.2) 

j^pd m — 1 z 

Note that as to — > 1, the first term in the integrand of converges to plogp which we refer to as 
the TO = 1 case. Using above energy structure, the existence and uniqueness properties of Eq. (|l.ip . 
in some appropriate Sobolev space, has been obtained in |BS| (also see [S] and [BRBj for relevant 
results). 

Compared to the well-posedness theory based on energy methods, few results has been known for 
pointwise behaviors of solutions, due to the lack of regularity estimates: the difficulty for regularity 
analysis lie mainly in the fact that the solutions are not necessarily positive (i.e., strictly positive) 
due to the degenerate diffusion. This is what we address in the first part of our paper. In addition, 
in the non-compact setting, the plausible limiting solutions tend to be trivial; here, since mass is 
conserved, even in the "worst" of cases, there is always the uniform stationary state. Most of the 
rest of this work is concerned with the approach to the asymptotic state. 

• Regularity properties Due to the degenerate diffusion, one cannot expect smooth solutions of 
Ea. (|l.l|) : even for (PME), Holder regularity is optimal, as verified by the self-similar (Barenblatt) 
solutions (see [V]). On the other hand, the solution of (PME) is indeed Holder continuous (again, 
see [V]), which motivates the question of Holder regularity of the solution of our problem Eq. (|l.ip . 

Note that, if we choose V as a. moUifier approximating the Dirac delta function, formally the 
nonlocal term approximates 

Therefore it is plausible that, at least when p is bounded from above, diffusion dominates when 
TO < 2 and the aggregation dominates when to > 2. Indeed we will show that, when to < 2, the 
effect of the aggregation term is weak enough that it is possible to locally approximate solutions 
of Eq. p.ip with those of (PME). As a result. Holder regularity of solutions of Eg. p.ip for to < 2 
follows. As for TO > 2, we show that solutions are continuous "uniformly in time", based on the 
result of Dibenedetto ( [Dibj ). For all to > 1, we also show that the L°° norm of solution is uniformly 
bounded from above depending on the and norm of the initial data (see Theorem 12. ip which 
is of independent interest. 

• Asymptotic behavior Our next result, partly an application of the first result, is on the asymptotic 
behavior of solutions of Eq. (|l.ip in the periodic domain T^. Wc work in a periodic domain because, 
primarily, wc are interested in finite volume problems and provides the most convenient boundary 
conditions. Even though asymptotic behavior for to < 2 has been studied before in various references 
(e.g., |S], and jHV| for a more singular interaction kernel) this is one of the first such result for these 
type of domains to the best of the authors' knowledge. One difficulty specific to the periodic setting 
is that the radial symmetry is not preserved over time, and thus exact (non-constant) solutions - 
always useful in these contexts - are not readily available. We also point out that in the case to > 2, 
there exist solutions which assumes zero value, possibly with compact support. Asymptotic behavior 
of such solutions are, in general, an interesting and difficult question, even for radial solutions in 
(see jKYj). 

Intuitively, one expects that when the diffusion term is "dominant" in Eq. (|l.ip . the solutions 
would converge to the constant solution as time goes to infinity. We show that this is indeed the 
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case when 1 < to < 2 and 9 is sufficiently small. However, with most interactions (specifically, V 
being not of positive type) there is a linear instability that sets in at some 0' ~ 0'^[ni) < oo which 
is determined by the minimal coefficient in the Fourier series of V (sec Section 4). It is not hard 
to show that for all to, when 6 > O'^, the functional in Ea. (|1.2p has non-constant minimizcrs (and 
the constant solution is not a minimizer - in fact, not even a local minimizer). However as has been 
shown explicitly for to = 1 under reasonable conditions - pertinently d> 2 - this "transition" occurs 
at some 9t < 9^ |CP| . Presumably, this argument holds in great generality. It is therefore somewhat 
surprising that for m = 2 the transition occurs exactly at 9 = 9K 

More precisely, for m 2, we show that for 9 < 9^ (the subcritical case), the constant solution 
is the only minimizer and is stable. Indeed we can actually show that for all bounded initial data 
p{x, 0), the dynamical p{-, t) will converge to the constant solution po exponentially fast in i^-norm. 
See Section 4 for detailed discussion on critical and supercritical case. When 1 < m < 2, the energy 
is no longer in the form of an L^-norm, and our Fourier-transform based approach does not generate 
a transitional value for 9. However, when 9 is sufficiently small, similar approach used by one of the 
authors in |CPj yields that the constant solution is the global minimizer. Moreover, we show when 
9 is sufficiently small, the solution uniformly and exponentially converges to the constant solution. 

Below we sketch an outline of our paper: In Section 2 we first give a uniform upper bound for 
the weak solution to porous medium equation with a drift for to > 1, then prove Holder continuity 
of the weak solution when 1 < m < 2. In Section 3 we apply the Holder continuity result to 
a nonlocal aggregation equation. In Section 4 we use Fourier transform approach to study the 
nonlocal aggregation equation when to = 2, and prove the exponential convergence of the weak 
solution in the subcritical case. Analogous results for 1 < to < 2 is established in Section 5. When 

1 < TO < 2, for 9 sufficcntly small, we prove there is also exponential convergence. 

2 Holder Continuity of the Solution of PME with a Drift 

In this section, we study the regularity of the porous medium equation with a drift, where the drift 
potential may depend on time: 

= A(p'") + V • (pV$) inl], (2.1) 

with Neumann boundary condition on dil. Here we may assume $7 is a bounded open set in R*^, 
where c? > 1, but all the results in this section certainly hold for periodic domain as well. We 
assume 1 < to < 2, the initial data p{x,Q) G L°°{VL) n L^{VL), the potential $(a;,t) S C{Q. x R+), 
and that $(-,i) G C^i^) for all t > 0. 

Before even stating the main result, we will first prove that p G L°°{il x M+). When $ does not 
depend on t, Bertch and Hilhorst in |BHj proved a uniform L°° bound of p by comparing p with 
an explicit supersolution which does not depend on t. When $ is a function of both x and t, using 
arguments similar to those in |KLj . wc aquirc an L°° bound for p which doesn't depend on t: 

Theorem 2.1. Suppose to > 1. Let p be the unique weak solution of Eq. l\2.1} with Neumann 
boundary condition, with initial data p{x,0) € L°° (Jl) Ci (il) . We assume that the potential ^{x,t) 
satisfies $(x,t) € C(f2 x M.~^), and <!'(•, i) S C^(J7) for all t with uniformly bounded norm. Then 
there exists M > 0, such that |l/o(-, i)||L=°(i2) 1^ for all t, where M depends on ||/o(2;, 0)||/^oo('f2), 
llp(a;,0)||ii(n), suptg[o^oo) ll*(-.i)llc2(o), andm. 

Proof. Wc begin with implementing the following scaling: Let 

p{x, t) = a™-i p[x, at), 
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where < a < 1. Let us choose a sufBcicntly small such that 



a<mm\{--——. )'"-\( „ , )""''^?^ 1' (2-2) 

where cq is a sufficiently small constant - certainly less than 1 - depending only on m and d and whose 
precise value will be determined later. By choosing a in this way, we have both |jp(a;, 0)||x,i(n) ^ 
Co < 1 and ||/5(a;, 0)||ioo(Q-) < 1 and, moreover, that p is a viscosity solution to the following PDE: 

Pt = A/5™ + V-(/5Vl>), (2.3) 

where $ := a<I>. From the definition of a we know ||$(-, t)||p2(f2-) < 1 for all t. 

Our preliminary goal is to show \\p{x, l)||L°°(n) ^ 1; then we can take p{x, 1) as the new initial 
data and iterate the argument to get a uniform bound for all time. 

We will introduce another variable v, which is bigger than p and is of order unity in x [0, 1]. 
Let V be the viscosity solution to the following equation 

Vt^V- {mv"'-^\7v + wV$), (2.4) 

with initial data v{x, 0) = p{x, 0) + . Since v solves the same equation as p with bigger initial 
data, we can apply the comparison principle for the porous medium equation with drift, which was 
established in Theorem 2.21 of |KL| . This comparison principle immediately implies v{x, t) > p{x, t) 
for all (x,t), hence it suffices to show |jw(-, l)||L=c(f2) < 1. 

One can check easily that v{x, t) [||w(-, 0)||Loo([2)]e-^* - where K := sup^^jQ q^-, ||<1>(-, i)||c2(n) - is 
a classical supersolution to Eq. (|2.4p and hence also a viscosity supersolution. Noting that the initial 
data of V satisfies, for all x, < u(x,0) < 1 + \e~^ ^ the comparison principle gives the following 

upper bound for v: 

\H-ML-in) < [|l«(-,0)|U^(o)]e^* < (1 + ie-i)e*. 
Similarly we can find a classical subsolution which gives the lower bound 

\\vi;t)\\L^m > [||«(-,0)|U»(o)]e-^' > ie-^e-*. 

Combining the two inequalities above, we have 

v{x,t) £ [ie"^,e + for allx ei},t e [0, 1]. 

We would like to refine the estimate above and get a better estimate at t = 1. By treating the 
diffusion coefficients mv"^~^ in Ea. (|2.4[) as an a priori function, - which we denote by b{x, t) - then 
we may say that v solves a linear equation of divergence form, where the diffusion coefficient is of 
(the order of) size unity: 

Vt^V ■ {b{x,t)Vv + vV^), (2.5) 

where h{x,t) := mv"''-^{x,t) G [m(ie-2)™-i, m(e + i)"-i] for a\\ x e fl, t e [0,1]. 

In particular, since Eg. p.Sp is linear, we can decompose v as vi + V2, such that vi solves Eg. p.Sp 
with initial data vi{x,0) = p{x,0), and V2 solves Ea. (|2.5p with initial data V2{x,0) = ^e~^. We 
claim that Vi{x, 1) and V2{x, 1) are both bounded by i, for all x & il. 

For vi, first note that due to the divergence form of Ea. (|2.5p . the norm of vi is conserved, 
i.e. ||wi(-, 1)||li(J2) = Cq. Since b is bounded above and below away from zero, then by |LSUj (see 
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Theorem 10.1, pp. 204), 1) is Holder continuous, where the Holder exponent and coefScient do 
not depend on cq, as long as cq < 1. So if we choose cq to be sufficiently small, we have vi{x,l) < i 
for all X G n. 

For V2, we can directly evaluate the necessary L°° bounds: 

supw2(a;, 1) < ell'^*ll°° supu2(x,0) < e l-e^^ = ^ 

(where again, on the basis of continuity, we may now talk about the supremum). 

Combining the two estimates together, we have sw^i^vix^ 1) < 1. which implies supj./3(a;, 1) < 1 
from our discussion above. Also, for < t < 1 we have p{x,t) < v{x,t) < e + 1/2. Then by treating 
p{x, 1) as initial data and iterating the same argument, we get sup^. p{x, t) < e + 1/2 for all t, i.e., 

p{x,t) < (e + ^)a"^ for all x eQ,t>0. 
Now plugging in the definition of a in the above and the bound becomes 

p(^,i)<(e + i)max{||p(x,0)|Uoo(o), ||$||^^)} inOx[0,oo). 

□ 

Remark 2.2. In the statement of Theorem 12.11 we assumed that 17 is a bounded open set, with 
Neumann boundary conditions. The same proof also applies to Dirichlet boundary condition. In- 
deed, the L°° bound we obtained is independent with the size of ft, and the same proof works as 
well when ft = 'R'^. However, ostensibly, the L°° norm of p should be of the order L"'' and, even if 
true in the initial data, we cannot establish that this order is preserved at later times. 

Since p{x,t) is uniformly bounded for all {x,t), DiBenedetto has shown in [Dib] that p{-,t) is 
continuous uniformly in t: 

Theorem 2.3 ( |Dib| ). For any m > 1, let p be the weak solution to jjJg. p.ip with initial data 
pix,0) e L°°{n) n L^{n). Let the potential $(a;,i) satisfy $(x,t) € C{n x R), e C'^{n) for 

all t, moreover supj ||$(-, (72(52 ■)] < 00. Then for all r > 0, p{x,t) is uniformly continuous in 
X [t, 00), and the continuity is uniform in x and t. 

Now we want to show when 1 < m < 2, for all t > 0, p{x,t) is uniformly Holder continuous in 
space and time in x [t, 00). Our main theorem of this section is stated as following: 

Theorem 2.4. Let 1 < m < 2. Let p be a viscosity solution of Eg. ^2.1^ . with initial data p{x,0). 
We make the following assumptions on p(-,0) and 

L ||p(-,0)||oc. < Ml and J^^p{x,0)dx < Mi. 

2. e C{n X R), and ||$(-, t)||c2(o) < M2 for some M2 > for all t > 0. 

Then for any < r < 00, u is Holder continuous m 17 x [r, 00), where the Holder exponent and 
coefficient depends on T,m,d, Mi and M2. 

Proof. To prove the Holder continuity of p, our goal is to show that for any (a;o,to) S x [r, c»), 

OSCB{xo,a^)x[to.to+a^]P < Ca'' (2.6) 
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for some C, 7 > not depending on a, for a satisfying < a < mini : v^j (where c is a 

2c 

constant to be determined soon). 

Bearing in mind that wc want to zoom in on the profile and look at the oscillation in a small 
neighborhood, it makes sense to start with a parabolic scaling with scaling factor a. Let 

p{x, t) := p{ax, a^t + (to - a^))), (2.7) 

and our goal Eq. (j2.6p would transform into 

< Ca\ (2.8) 

Here p{x, t) is defined in the domain Vt x [0, 00), where {x e R'' : ax e il}. and, it is noted, the 
early portion of the time domain had been omitted. We readily see that p is the viscosity solution 
to 

pt Ap" + V • (pV$) in X [0, 00). (2.9) 

Here, the initial data reads p(a;, 0) = p{ax,to — a^), which has an a priori L°° bound depending 
on m, d. Ml, M2 due to Theorem 12. II Moreover, in the above $(x, t) := $(aa;, a^t + {to — a^))) and 
hence |V$| is bounded by aM2. We wish to compare p with w, where w is the viscosity solution to 
the porous medium equation 

wt=Aw" in fix [0,00), (2.10) 

with initial data w(-,0) = p(-,0). Since Ea. (P3)) and Eq. (P?TU)) only differ by the term V • (/5V$), 
we would expect 

|p- w| < Ca-^ in O X [1,2], (2.11) 

for some C>0,0</3<1 depending on m, d, A/i, Af2- 

The main part of this proof will be devoted to proving Ea. (|2.11[ ) is indeed true. Without loss 
of generality, we can assume that p{x, t) is a classical solution. First, if the initial data p{x, 0) is 
uniformly positive, then p{x, t) will be a classical solution for all time. This is because p will stay 
positive for any time period [0, T] (since inffj^ jg p{x, t) > exp{—t supjgjg.T] II ^^11 00) inf^g^ p{x, 0)), 
which implies that Ea. (|2.9p is uniformly parabolic for t £ [0,T] and hence the weak solution p is 
classical. 

For general initial data p(a;, 0), we can use approximation as follows. Let p„ and w„ solve Ea. (|2.9p 
and Ea. (j2.1ip respectively with initial data p{x,0) + 2^"; n sufficiently large. As discussed above, 
Pn would be a sequence of classical solutions. If we can obtain |/5„ — Wn\ < Ca^ for all n, (where C, /? 
doesn't depend on n), then Ea. (|2.1ip would hold for p and w as well, since as n — >■ 00, comparison 
principle yields pn{x,t) \ p and Wn(x,t) \ w uniformly in x,<. 

Note that one cannot directly compare p with w, due to the fact that the term V • (/5V$) contains 
Vp • V$ and hence docs not have any a priori bound. In order to bound this term, it will help to 
change from the density variable p to the pressure variable u. Let 



TO — 1 

then Eq. (|2.9p becomes 

iLt^im- l)uAu + |Vwp + Vm • V$ + (to - l)uA$, (2.12) 

which will enable us to use jV-up plus a constant to control the term Vu • V$: Recall that |V$| < 
aM2, which gives us the following bound 

|Vu • V$| < aAfalVwl < a[|Vw|^ + ^(Ms)^]. 
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Also, due to the fact that (m — l)u{x,t) < Ci in f2 x [0,2], (where Ci, which depends on m,d,AIi 
and M2, is related to the L°° bounds on p) we obtain 

|(m - l)uA$| < a^CiAfz < aCiA'h- 

Putting the above two bounds together, and by choosing c such that c > C1M2 + (il/2/2)^, u will 
satisfy the following inequality 

ut > (m-l)uA{t+(l-ca)|V?ip-ca for all x e n,t e [0,2]. (2.13) 

Note that we assumed a < (2 — m)/{2c) in the beginning of the proof, we have ca < (2 — m)/2. 

In order to make Eq. (|2.13p look similar to the porous medium equation in the pressure form, we 
apply the rescaling ui = (1 — ca)u. Then ui satisfies 

(wi)t > (m" - l)wiAwi + |Vwi|^ - ca(l - ca) for all a; e f^, t e [0, 2], (2.14) 

where 

Tn — 1 

:= h 1. (hence ca < (2 — m)/2 implies that 1 < < 2) (2-15) 

1 — ca 

Now Eq. ((2?T4)) has the same form as the porous medium equation in the pressure form, minus an 
extra constant term ca(l — ca). To take advantage of the existence and regularity results for equations 
with divergence form, we change the pressure variable back to the density variable (however here 
the power is m~ instead of m), i.e., we define pi such that 



or in other words, 



(1 — ca)u = ui = — p^ ^, (2-16) 

m — 1 

PI = {—r-p'-"^ = (-i±^)""pi-. (2.17) 

m 1 + ca/m 

Due to the positivity of u, we know pi is positive as well. Hence when we plug Ea. (|2.16p into 
Ea. (|2.14p . after canceling a positive power of pi on both sides, we obtain 

{pi)t > Ap'i"" - ca(l - ca)pl-"'' in n x [0, 2], (2.18) 

Note that the term ca(l — ca)p\^^^ has an a priori upper bound: since 2 — m~ > and pi is given 
by Eq. (|2.17p . we have c(l — ca)p^~™ < M, for some constant M depending on to, d, Mi, M2. 
Let us denote by p~ the weak solution of 

{p-)t = A{p-\p-r--')^Ma, (2.19) 
with initial data the same as pi{x, 0), which is 

p-(a;,0) = ( — )""p(z,0)i— (2.20) 



Since f2 is a bounded domain, we have Ala G L^(J7) for all p > 1, and the existence of weak 
solution of Ea. (j2.19p is guaranteed by Theorem 5.7 in jV|. That theorem also gives us a comparison 
result that, a.e., pi > p^ ■ 
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Moreover, note that the "a.e." above can in fact be removed, since both p and p" are continuous 
in f2 X [0,2]: the continuity of p is given by Theorem 12. 3[ and the continuity of p~ is given by 
Theorem 11.2 of |DGVj . Therefore we have the following comparison between p~ and p: 



< (^)" p'"™ in n X [0,2] (2.21) 



m 



Since m/m~ = 1 + 0{a), and p is bounded in x [0, 2], Eq. (|2.21|) implies that p — p^ > — Cio, 
where Ci depend on m, d, Mi,M2- 

Analogous to the definition to p~ , we define p+ to be the weak solution of 

(p+)t = A((p+r+) + Ma, (2.22) 

with initial data 



m 



p+(x,0) = (—)■" p(x,0)^+™, (2.23) 
where 

fn — 1 

m+ := h 1. (hence 1 < ?7i+ < 2) (2.24) 

1 + ca 

Then analogous argument would lead to p — p^ < Cia. Summarizing, we have obtained 

p- -Cia<p<p+ + Cia in fie [0,2], (2.25) 

where Ci depends on m, d, Mi, M2- 

To prove Ea. (|2.11l) . it suffices to show \p^ — w\ < 0{a^) for some /3 > 0, which is proved in the 
following lemma. 

Lemma 2.5. Let 1 < m < 2. Let w be the viscosity solution of the porous medium equation 

wt = Aw" inVty. [0, oo) (2.26) 

where the initial data w(x, 0) satisfies i(j(x,0) — p(x, 0). 

Let p~ and p'^ be the weak solutions to i?o. (|2.19p and i?q. (|2.22[) respectively, where < a < 
(2 — m)/(2c) is a small constant, and the initial data is given by i?o. (|2.20p and i?o. (|2.23p . T/ien 

- wl < Ca'' in fix [1,2], (2.27) 

where C and /S depends on d, to, Afi, M2. 

The proof of Lemma 12.51 is the content of the appendix in Section [S] Putting Lemma 12.51 and 
Eg. ((2^ together, we obtain Eq. ((2JT|) . and we will use this to (immediately) prove Eq. (j2.8p . 

Since w solves the porous medium equation. Theorem 7.17 in |V] gives us the Holder continuity 
of w: 

OSCBix,a)xli,i+a^]W < Ca" , for aU x e Cl, (2.28) 

where C and a depends on ||?ii(-,0)||oo (and hence depends on m,d, Mi, M2). 
By putting Eq. (|2.29p and Eq. (|2.1ip together, we obtain 

OSCB{x,a)xli,i+a'^]P < Co'' , for aU X e n, (2.29) 

where C depends on to, d. Mi, M2, and 7 — minjo;, /?} (hence also depends on m, d, Mi, M2). Hence 
Eq.^^ is proved. □ 
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Remark 2.6. For m > 2, Holder continuity of the solution to Ea. (j2.ip is still open. Indeed, 
concerning the present approach - which closely parallels that of [KL| . [K] - when m > 2 we have 
that m~ = l + (m— 1)/(1 — ca) > 2. Hence the "inhomogeneous" term in Eq. (|2.18p . which is 
proportional to pt^"™ ), would actually be divergent in places where p — > 0. This indicates that 
another approach will be required. 

3 Application to Aggregation Equation with Degenerate Dif- 
fusion 

In the following two sections, we study Eg. p.ip in the domain T^, the d-dimcnsion torus of scale 
L. Here is a non-negative constant, and, of course, * denotes convolution in T^. We make the 
following assumptions on V{x): 

(VI) V{x) = V{-x) for aU xeTf. 

(V2) V{x) e C2(T^), with \\V{x)\\c2f_jd-j = C for some constant C < oo. 

Moreover, we have in mind : M'' R compactly supported with the diameter of the support 
smaller than L. In particular we do not envision "wrapping" effects and Jj.d \V\dx may be regarded 
as independent of L. 

Our goal in this section is to show the Holder continuity of the weak solution to Ea. (|l.ip for 
1 < m < 2, and uniform continuity of the weak solution when m = 2. First, we state the definition 
of weak solution to Ea. (|l.ip and a existence theorem from |BSj . 

Definition 3.1 (Weak Solution). Let m > 1, and let us assume that p{x,0) is non-negative, with 
p{x,0) G L°°{T'j^) and consider a potential V that satisfies the assumptions (VI) and (V2). 
A function p : Tf x [0,r] [0, oo) is a weak solution to Eq.l^^ if P ^ L°°(T^ x [0,r]), 
e L^{0,T,H\T[)) {i.e., ||/9(-, OIIhmt^ ^ i^(0,T)) and pt G L^{0,T, H-\Tl)) and for all 
test function (p £ H^{T'l), for almost all t £ [0,T], 

< pt(t), (j)>+ [ Vip'^it)) ■ V0 + 0L''(2-™)p(t)(VF * p{t)) ■ V<pdx = 0. (3.1) 

In |BS| , existence and uniqueness of weak solution are proved: 

Theorem 3.2 (Bertozzi-Slepcev). Let m> 1 and consider V that satisfies the assumptions (VI) 
and (V2). Let p{x,0) be a nonnegative function in L°° (Tf^) . Then the vroblem i?o. (jl.ip has a unique 
weak solution on x [0,T] for all T > 0, and furthermore p £ C(0, T, (T^)) for all p £ [l,oo). 

By treating * 1^ as an a priori potential, we can apply our results in Section 2, and 

obtain L°° bound of p which does not depend on T, together with uniform continuity of p, and 
Holder continuity of p for 1 < m < 2. 

Theorem 3.3. Let m > 1 and consider V that satisfies the assumptions (VI) and (V2). Let 
p{x,t) he the unique weak solution to Ea. l\l.l\ given by Theorem \3.2l with nonnegative initial data 
p{x,0) £ C(T^), which satisfies J^.^ p{x,0)dx = 1. Then ||/o(-'J;,i)||L°°(T£x[o oo)) bounded, where 
the bound only depend on sup^p(x,0), 9, \\V\\c2 and L. 
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Proof. To begin with, note that Theorem 13.21 guarantees the existence and uniqueness of the weak 
solution to Eq. (|l.ip . which we denote by by p. Now we treat $ := 6IL''(2-™)p ^ y as an a prion 
potential, and we obtain the following estimate of $ assumption (V2): 



m-McHTi) < OL''^'-^-^\pi;t)\\LHTi)\\V\\cHTt) 

We denote by pi the unique weak solution to the equation 

= Apr + V-(piV$) (3.2) 

with initial data pi{-,0) = p{-,0), where the existence and uniqueness is proved in |BH] . Theorem 
12.11 implies sup^ bounded uniformly in t. Moreover, note that p also satisfies the weak 

equation for Ea. p.2p . hence p must coincide with pi, which yields a uniform bound of p which 
doesn't depend on time. □ 

Applying Theorem 12.31 to Eq. p.2[) . we have the continuity of p uniformly in t for m > 1 ~ in 
particular (in light of Theorem 13.51 below) for the case to = 2. 

Theorem 3.4. Let m > 1 and consider V that satisfies the assumptions (VI) and (V2). Let 

p{x,t) be the unique weak solution to £'o. (|l.ip given by Theorem \3.2l with nonnegative initial data 
/9(-,0) satisfying ||/o(-, O)||^=o(-Td) < oo, and ||/3(-, 0)||^i(Td-) = 1. Then for any t > 0, p is continuous 
in T'l X [r, oo), where the continuity is uniform in both x and t. 

Proof. Follows immediately from the above reasoning. Theorem 13. 31 and Theorem l2.3l □ 

Applying Theorem [2^ to Eq.(|3.2|). with $ 6L'^^'^^"^^p * V wc have the Holder continuity of p 
for 1 < TO < 2. 

Theorem 3.5. Let 1 < to < 2 and consider V that satisfies the assumptions (VI) and (V2). Let 

p{x,t) be the unique weak solution to Ea. l\\.l\ given by Theorem \3.2[ with nonnegative initial data 
p{x,0) satisfying \\pi-,0)\\L'^(jd-^ < oo, and \\p{-,0)\\l'-{T'1) = 1- Then for any t > 0, u is Holder 
continuous in x (r, oo), where the Holder exponent and coefficient depend on T,m,d,9, L and C 
and the L°° norm of the initial condition. 

Proof. Follows immediately from the preceding reasoning, Theorem 13.31 and Theorem 12.41 

□ 



4 The Case m = 2: Analysis Via Normal Modes 

In this section, we will use Fourier Transform to study the PDF in Fq. ljl.ip . and this method works 
best when to = 2. We continue to assume, without loss of generality that ||p(2:, O)!!^!^^^ ) = 1, 
however from the perspective of functional analysis, the homogeneity of the special case to = 2 
makes even this stipulation redundant. 

The dynamics in Eq. (|l.ip is governed by gradient flow for the "free energy" functional 

Mp)^ I p" + \ep{p*v)dx. (4.1) 
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For the analysis of the functional J-g, since we are assuming p(x,0) integrates to 1, we shall denote 
by ^ the class of probability densities on which also belong to L^(T^), i.e. 

^ := {/ e Li(Ti) n L^Ti) : = 1}. (4.2) 

Special to the case m = 2 is that the functional J-g{-) can be expressed in a simpler form if we 
express p in terms of its Fourier modes. We write 

p{k) = / p{x)e-'''''dx 
where k is of the form k — with n £ Z'^. With these conventions we have 



k 



and, in terms of these variables, Eq. (|4.ip becomes 

-^«(/') = idEl/5(fc)l'(l + ^^nfc))- (4.3) 

k 

On the basis of Ea. (|4.3p . a salient value of 9 emerges: We denote this value by 6'", which is defined 
via 

[0«]-i:=imax{|l>(fc)|;f(fc)<O}. (4.4) 

Formally 6*" may be designated as +oo in case V{k) > for all fe ^ - i.e. if V is (essentially) 
of positive type. For the purposes of the present discussion, we shall assume otherwise. Different 
values of 6 separate our problem into 3 cases: 

1. (subcritical) When < 6*", we have 1 + ^9V{k) > for all k e Z^, then under the restriction 
p(0) = 1, it is manifest that global minimizer for J-e{p) in ^ is the constant solution 

If 1 

Po{x) := j-j J^^ p{x,0)dx = j^. (4.5) 



2. (critical) When 6 = 6'", we have still have 1 + ^9V{k) > for all A: e however now 
there is a set Kl* (containing at least two elements) defined by the condition that for /c G K", 
1 + ^9'^V{k) = 0. In this case the global minimizers for J-e{p) in ^ take the form 

p{x)^po+Y.Cke^'-\ (4.6) 

fceK« 

where C-k = Ck and, of course, subject to the restriction that the resultant quantity is non- 
negative. 

3. (supercritical) When 9>9'^,we have 1 + ^9V{k) < for some k G Z'^. In this case the constant 
solution pq is not even a local minimizer of J^e in let alone global minimizer. 
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Remark 4.1. The above ~ which is manifest for m = 2 - is in sharp contrast to the cases m ^ 2. 
In particular, for general m there is an analogous quantity 6^ given by 



:=-max{|f(fc)|;f(fc)<0} 

where items (1) - (3) are suggested. However, the following was shown for m = 1 and, presumably 
holds for all m ^ 2: While for 6 < 6\ the constant solution has "some stability" (c.f. |CP| Theorem 
2.11 for the case to = 1) there is a 0^ < where global considerations come into play. In particular, 
aX 9 = O-T, there is a non-uniform minimizer for J-g^ (•) which is degenerate with the uniform solution. 
Moreover, for 9 > 9^ (which implies, in particular, at 9 = 9"^) the uniform solution is no longer a 
minimizer. 



4.1 The subcritical case, when m=2 

In the subcritical case, the constant solution pq is the only global minimizer of J-g in Our goal 
in this section is to show for every non-negative initial data p{x,Q) G L°°(T^) which integrates to 
1, the weak solution p{x,t) converges to pq exponentially in L'^{Tf) as < — >■ oo, where po is as given 
in Eq. (|4?5|) . 

By formally taking the time derivative of the free energy functional, a simple calculation indicates 
that e.g., at least for classical solutions to Eq. (|l.ip . the free energy is always non-increasing: 

^Fe{p) = -I p|v(^^p™-i+0L'^(2-™V*^)|'rf^- (4.7) 
at Jjd^ TO — 1 

In [BS| . it is proved that Ea. (|4.7p is indeed true in the integral sense: 

Lemma 4.2 (Bertozzi-Slepcev). Consider V that satisfies the assumptions (VI) and (V2). Let 

p{x,t) be a weak solution of in x [0,T]. Then for almost all r e [0,r], 

Mpi; 0)) - Mpi;r)) > r f p|V(^^p'"-i + eL-^^^-rn)^ , V)\'dxdt (4.8) 

Jo Jri - 1 

Remark 4.3. Theorem 13.41 implies that p{-,t) is a continuous function of t, hence J-g{p{-,t)) is 
continuous in t as well. Therefore Ea. (j4.8|) indeed holds for all r G [0,T] and, moreover, Ea. (|4.7p 
may be regarded as a differential inequality. 

In the following lemma, we show when 9 < 9'^, the free energy will decay to the free energy of the 
global minimizer as t — > oo. 

Lemma 4.4. Suppose to = 2 and consider V that satisfies the assumptions (VI) and (V2). 
Further suppose that 9 < 9^, where 9^ is as given in j?q. (j4.4p - including 9^ = oo if V is of positive 
type. Let p{x,t) be the weak solution to i?a. (|l.ip on [0, oo) x T^, with non-negative initial data 
p{x,Q) E L°°(T^) which integrates to 1. Then J^e{p) ^ J^e{po) as t oo, where po is the uniform 
solution (as given in 

Proof. By Lemma we know Tg{p{t)) is a continuous and decreasing function of t, whose limit is 
bounded below by J^g{po), since po is the global minimizer of J^g in ^ when 9 < 9K Hence we can 
send r to infinity in Ea. (|4.8p . which gives 

[ [ p\V{2p + 9p *V)\^dxdt < oo. (4.9) 
Jo JtI 
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Then there exists an increasing sequence of time (tn)5^i, where hni„^oo tn — oo, such that 

hm / p{x,tn)\y{2p{x,t^) + ep{x,t,,)*V)\^dx = 0. (4.10) 



To avoid clutter, in what follows, we shall abbreviate p{-,tn) by p„. Recall that Theorem 1 2 . 1 1 gives 
us a uniform bound of ||pn||L°°(R<i)- In addition, by |Dib| . (p„) is uniformly equicontinuous, hence 
Arzela-Ascoli Theorem enables us to find a subsequence of p„ (which we again denote by p„ for 
notational simplicity), and a continuous function poa, such that 

lim |1p„ - poo|lL-(Tf ) = 0, (4.11) 
Wc next claim that || Vpf/^||^2(-xd ^ is bounded uniformly in n. To prove the claim, we first note that 

/ \^^p''J^+pl/^^iOpn*V)fdx= [ Pn\2Vpn+V{epn^V)\^dx^0. (4.12) 



To obtain the uniform bound for V pf/^ , due to the triangle inequality, it suffices to prove a uniform 
bound for pl/^V{9pn * V), which is true since pn is uniformly bounded in n and ||T^||(^2(t'' ) < oo 
due to (V2), hence the claim is proved. 

As a consequence of the claim, we obtain weak convergence of Vpf/^ in (along another sub- 
sequence) And, it is clear, the limit is just Vp^^ due to the uniform convergence of the (pn)- 
(Moreover, this places Vpi^^ £ L^{Ti : W^)). Thus: 

^ Vp^2 as n oo weakly in L^{Ti : M"^). (4.13) 

Let 

Bn-.^^vpl^' + pU^yiepn^v). 

Then Ea. (|4.1ip and Eq. (|4.13p and an additional uniform convergence argument identifying the weak 
limit of pl/'^V{9pn * V), implies that Bn weakly converges to Boo in L^, where 

Boo := ^VpL^2 ^ pi/2v(0^^ , V). 

On the other hand, recall that Ea. (|4.12p gives us that i3„ — )■ strongly in L^, thus we have Bqo is 
indeed i.e.. 



/ \^^pT + Po^''^iSPoo*V)fdx^ [ poo\2Vpoo + Viepoo*V)\^dx = 0. 
Jjd 6 Jfd 



(4.14) 



In particular, then, V(poo + 56*^00 * V) is zero a.e. on the support of poo- Now p^o certainly admits 
a weak derivative which, clearly, is non-zero only on the support of poo- Thus, from the preceding, 
we can write 

/ Vpoo • V(poo + ^rfilpoo *^^)dx = 0. (4.15) 

Now, we wish to express the above as a Fourier sum which requires some additional justification. 
To this end we claim that p^o is Lipschitz continuous - i.e., in Vl^^'°°(T^) - which places both entities 
in L^(T^) and vindicates the use of explicit formulas. 
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The equation Vpoo = — ^0V(y * poo) valid on the support of p shows that in the various com- 
ponents where poo is positive, it is at least C^. Indeed, in general. Hypothesis (V2) immediately 
implies ||poo(a;) * ^IIc2(t|^) — IIPooUli II^IIc2(t^) so whenever poo satisfies this (m = 2 version of the 
Kirkwood-Monroe) equation, we have Lipschitz eontinuity with uniform constant. We shall denote 
this constant by k. Now suppose that x,y & have Poo{x) and Pociu) positive. Let us assume, 
ostensibly, that x and y belong to different components. On the (shortest) line joining x and y, let 
Zx denote the first point, starting from x that is encountered on the boundary of the component of 
X and similarly for Zy. Then 

|Poo(a;) - Pooiy)\ = \Poo{x) - Poo{Zx) + Poo{Zy) - Poo{y)\ 

< \Pooix) - Pooiz^)\ + \PooiZy) - Poo(y)| 

< k[\x - z^\ + \y - Zy\] < k\x - y\; (4.16) 

the first inequality due to Poo{zx) = Poo{zy) = and the last inequality because all four points lie in 
order on the same line. A similar argument can be used if, e.g., Poo{x) is positive and Poo{y) is zero. 
All of this establishes enough regularity to unabashedly express Eq. (|4.15p in Fourier modes: 

= E^I'5oo(fc)P(l + ^^V^(fc))- (4.17) 

k 

By the defining property of 6^ we have l + }^9V{k) > for all fc 7^ 0, thus Ea. (|4.17p implies poo{k) = 
for all k ^ 0, i.e. Poo = Pa- 

Now, we may use the monotonicity in time of Tg{p{t)) and we finally have 

lim Fg{p{t)) = lim Fe{pn) = J^eiPoo) = J^eipo) 

which is the stated claim. □ 

By combining the above result with the uniform continuity in time, we can show the solution will 
become uniformly positive after a sufficiently large time. 



Corollary 4.5. Under the assumption of Lemma \4.4\ we have 

lim \\p{-,t) - PoIIloo(t^) = 0, 



hence there exists T > depending on 9, \\V\\(j2(jd-^ and p(-,0), such that p{x,t) > po/2 /or all 
X £ Ti,t> T. 

Proof. We prove the statement in the display. Supposing that this is not the ease. Then there is a 
sequence of times, (t„) and points [yn) - yn £ - and a. S > such that 

\p{yn,Tn) - Pol > S. 

Now, going to a further subsequence, we have ?/„ — > yoo (with y^c G TT^ by compactness). But, along 
yet a further subsequence, not relabeled, we have, according to the arguments of Lemma 14.41 that 
p(-,T„) is converging uniformly and the limit must be po- Thus 

lim p„(y„,T„)= lim [p„(y„,T„) - p„ (2/00, T„)] + lim p„(yoo, t„) = po 
in contradiction with the preceding display. □ 
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Theorem 4.6. Suppose m = 2 and 9 < 6^ , where 0^ is as given in £'o. (j4.4p . Consider V that satisfies 
the assumptions (VI) and (V2). Let p{x,t) be the weak solution to i?q. (|l.ll) on [0,oo) x T^, with 
non-negative initial data p{x, 0) € L°°(T^) which integrates to 1. Then J-g{p{t)) decays exponentially 
to J^s{pq), where the rate depend on p{x,Q). Moreover, \\p{-,t) — po||L2(Tf) ^ exponentially, i.e. 

< Mp{t)) - Mp^) < Ci exp(~^i), 

and ^ 

\\p{t) - Po|1l2(tJ) < Ca exp(-^i), 

where d and C'l and C2 depend on 6, V and p(-,0). 

Proof. By Lemma l475l there exist some T > depending on 6, V and 0), such that p{x, t) > po/2 
for all X eT'l,t > T. Then for all t2 > ti > T, Eq.dH)) becomes 

Mp{;ti))-Mpi;t2)) > f ^\^{2p + ep*V)\^dxdt 

t: 



pt2 

= Poc {Tg{p{-,t))-MPo))dt, (4.18) 
Jti 

where c' = 2minfc^o + ^^^^s)): which is positive when 9 < OK 

In the spirit of Remark 14.31 we may regard the above as a differential inequality for g(t) := 
T0{p{-,t)) — Tg{po); the inequality reads 

-^>Poc'g{t). 

This immediately integrates to yield g{t) < g{T)cxp{— poc' {t — T)} for t >T. I.e., 

J^{p{-,t))-J^{po)<Cc-P°-'K 

Since F0{p{-,t)) — Te{po) is comparable to \\p{t) — polU^, we have \\p{t) — po\\l2 — > exponentially 
with the same rate. □ 

Remark 4.7. It is remarked that, via comparison to linearized theory, the above is essentially 
optimal. (The results differ by a factor of two which comes from the definition of T =: Ti/2. Using 
Te = sup{t > I ||p(-,^) — Po||L°°(Tf) > £Po}i the long time asymptotic rates are actually in complete 
agreement.) Moreover, while for L of order unity, the result stands: c' - with or without an additional 
factor of two - might well be optimized at a wave number of order unity. However, as L — >■ 00, it is 
clear that 

min|fcp(l + if(fc))^(^)2(l + if(0)). 

So, in particular, for large L the rate scales as L^('^+2) - a result which may be an artifact of our 
normalization. 



15 



4.2 Some remarks on the supercritical case, when m = 2 



When e > e\ we have 1 + \eV{kQ) < for some kQ — ^fiQ, where fiQ £ Z"^. In other words, at least 
one of the coefficients of the free energy Ea. (|4.3p is negative. In the next proposition we show that 
in this case the constant solution po is not linearly stable. 



Proposition 4.8. Suppose m = 2 and < 0* , where 6* is as given in Consider an 

interaction V that satisfies the assumptions (VI) and (V2). Then the constant solution po is not 
a local minimizer of the free energy j?g. (|4.3|) in 



Proof. We choose fco — such that l + ^9V{ko) < 0, where no e Z^. Wc add a small pertubation 
erj to the constant solution po, where 



Then 



27rno • X , 
V ■■= cos( ). 



which is strictly less than J-e{po) by the defining property of fco . □ 

Remark 4.9. In fact, using the same perturbation term in the proof, we would know that when 
9 > 9'^, any strictly positive function is not a local minimizer of the free energy Eq. (|4.3p . 



In the supercritical case, while Eq. (|4.3p immediately implies that po is not a local minimizer of J^e 
in it gives us little information about what is the global minimizer. The difficulty comes from 
the restriction p(x) > for all x, which evidently plays an important role in the supercritical case, 
since any minimizer should touch zero somewhere due to Remark 14.91 After Fourier transform, the 
non-negativity of p actually gives us infinite numbers of restrictions, which causes the difficulty. 



5 Exponential decay for 1 < m < 2 and weak interaction 

In this section, we continue our study of Eq. ljl.ip with m e (1,2) and here we will assume that 9 
is "small" . Unfortunately, 9 will not be uniformly small in volume. In particular, we shall require 
f)ljd.{2-m) small number of order unity and, under these conditions we shall acquire all the 

results of the previous section. We claim that without additional (physics based) assumptions - in 
particular H-stability of the interaction - the above condition is essentially optimal. Specifically, 
our cornerstone result of a unique stationary state does not hold for non-H-stable interactions when 
is a sufficiently large number of order unity. However, from an esthetic perspective, this 
uniqueness result is the sole instance where 9L'^^'^~"^'> must be considered small. In the aftermath of 
Proposition 15.11 and its corollary, we will only require 9 itself to be a small quantity. 

We start with a priliminary result (which is, actually, just a quantitative version of the argument 
used in Lemma [4.41 in the vicinity of Eq. (j4.16|) ). 

Proposition 5.1. Consider an interaction V that satisfies the assumptions (VI) and (V2). Let 



be a sufficiently small number of order unity. Let p denote any solution to the Kirkwood-Monroe 
equations which here read, whenever p > 0, 

m 
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and let 

R ■= ll/'llL~(Ti)- 

Then if Sq is a small number of order unity then R is also a small number of order unity (if L is 
large). In particular, 

R < K4 max{[£o] k^-d+i , L~'^} 

with K4 a constant of order unity. 
Proof. From the mean-field equations, 

iVp^-il < !^£o / |VF(x - yMy)dy < !I^\\V\\c,eo ^iSo. 
m Jjd m 

Let xo mark the spot where p achieves R. Then, for all x, 

> - Kieoix - a;o|. 

Thus, if r is the length scale of the region about xq where exceeds, a.e., ^i?™"^ we have 

r > 

2Kieo 

provided the right hand side does not exceed L. Otherwise, obviously, r ^ L. Since p integrates to 
unity we have, assuming r < L, 



^(i(m-l) + l ^ 2 



, pdx > K2r'^R > K2- -r -r-;R'^'^"'-^^+^ 



(with K2 a geometric constant of order unity) and otherwise we acquire the mundane bound. After 
a small step, the stated bound is obtained with an appropriate definition of K4. □ 

With the above in hand, we can establish that po is the unique stationary solution. We start with 



Corollary 5.2. Under the conditions stated in Provosition 15.11 if Eq is sufficiently small - but of 
order unity independent of L - the unique solution to the mean-field equations is p = pQ. 

Proof. From the mean-field equation, we may write 

Jjd m 

By recapitulating the Lipchitz continuity that was featured in the vicinity of Eq. (j4.16p we have full 
justification to manipulate classically under the integral. Letting iJ^^ denote the upper bound on 
the norm of p that was featured in Proposition l5.ll Then, pointwise a.e. on the support of p. 



since, we remind the reader, 2 — m > 0. In other words, 

0> / -l^|Vpp + -Vp-V(p*F)dx. 

J-Td RL m 
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Wc can again go to Fourier modes and the above reads 

o>^fc^ip(fc)n-^ + ^nfc)]. 

For Eq sufficiently small (but of order unity independent of L) the coefficient of is positive 

for all terms so the later must vanish identically. The desired result is proved. □ 

Based on the fact that po is the unique stationary solution, in the next lemma we prove that t) 
will converge to po uniformly, (but not with a quantitative estimate on the rate.) 

Lemma 5.3. Suppose the conclusions in CoroUaru \5.3\ are satisfied. Let p{x,t) be the weak solution 
to Eg. 1^1.1^ on [0,oo) x T^, with non-negative initial data p{x,0) G L°°(T^) which integrates to 1. 
Then sup^. t) — pq\ ^ as t ~^ oo. 



Proof. This is more or less identical to the proof of Corollary 14.51 based on Lemma 14.41 so we shall 
be succinct. Assuming the result false, we could find a sequence of times t„ — ?> oo and points 
Xn — > S such that p(-,i„) converges uniformly and yet |/o(a;„,T„) — pol > 5. So, denoting by 
Poo(') the uniform limit, we would have |poo(a;oo) — Pol > 

Hence, since pao is continuous, it is definitively not equal to po- However, any subsequential limit 
must satisfy the mean-field equation and by Corollary 15.21 this is uniquely po in contradiction with 
the preceding. This completes the proof. □ 

In the next lemma, we show that once p and po becomes comparable, J^eip) — J^eifo) also becomes 
comparable with L''^^~™^||p— Po||l2(t^)- Indeed, as alluded to earlier, this will be proved under the 
weaker assumption that 9 - not „ jg gj^all. We start with: 

Lemma 5.4. Suppose that 6 > is sufficiently small (but of order unity independent of L). Let p 
be such that ||p — poUt^ < ^Po- Then we have 

aL^(2-™)||p-po||i.(T.) <.F,(p)- J-e(po) </?i'^(^-'"'||p-po|li.(T.^ (5.1) 
for some a, /3 > of order unity. 
Proof. First, by any number of methods we have 



p{x)p{y)V{x - y)dxdy > -AVilp - Po|li2(Ti ); 



e.g., we may take, using the Fourier decomposition, Ky = l^"]^^- Similarly for a corresponding 
upper bound with a positive constant. Let us turn to the entropic-like terms. 

Writing p = po(l + v)i our assumption implies that I77I < i. From this it is easy to verify that, 
pointwise, 

(1 + r;)™ > 1 + mr; + (^)^-'"r;^ := 1 + mr; + mf , 

and for the other direction we have 

(1 +??)"< 1 + m7? + ""^"^ ^ ^) (3)^""?/^ := 1 + 77177 + 6772. 
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Thence p'" - + 77)'" - 1] = + 77)'" - 1 - 77777 + mil] > pj," [77777 + 077^]. So 

and similarly we have 

Combining this with the bounds on the energy term, the stated claim has been estabhshed. 



□ 



Finahy, in the next theorem, we prove that the free energy decays exponentially to its minimum 
value. 



Theorem 5.5. Suppose the conclusions acquired in Corollary [57^ are satisfied and suppose that 9 is 
a sufficiently small number which is of the order of unity. Let p{x, t) be the weak solution to Eq. Hl.l^i 
on [0, 00) X T^, with non-negative initial data p{x, 0) € L°°(T^) which integrates to 1. Then To{p{t)) 
decays exponentially to J-g{po). More precisely, 



J-(p(-,i))-J-(po)<Cie-^"" 



(5.2) 



for various constants c' and Ci. Similarly for the L^-norm of (p — pq) with a different prefactor. 

Proof. According to Lemma [5731 there exist some T > depending on 9, L, V and p(-, 0), such that 
\p{x, i) — Pol < ^Po for all a; G T^, t > T. Then for all <2 > > T, we manipulate the integrand on 
the right hand side of Ea. (|4.8p - the lower bound on Te{p{-,ti)) — J'e(p(-, ^2)): 



p|V 



777—1 



2'^ 777 - 1 



P 



-1|2 _ |6lid(2-m)Vp^^|2 



2 ^ 



> 



5Pr-'|Vpp - ^Poe^L2^(^-™)|V(p* F)p 



dx 
dx 
dx 



(5.3) 



where the value of g - which is always of order unity - depends on whether 2777 — 3 is positive or 
not. Note that aU terms are proportional to Pg™^^ = p^^^L'^'-^^"^\ 

Going to Fourier modes, the final (spatial) integral in the above string becomes 



Po 



rn — lj^d{2 — ni) 



where, for sufficiently small 9, we may assert that the summand is positive. 
We thus have 



Fe{p{-M))-^e{p{-M)) > 



ni—1 n t 

Po Pc 



U'^'''-\p- pofdxdt > p'^-^' [Fe{p{-,t))^Fe{po)]dt (5.4) 
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where in the above, (3 is the constant from Lemma 15.41 which has been conveniently absorbed into 
the definition of c': 

c'P ■.= m^n[k\g~lf\V{m] 

and in the final step we have used Lemma 15.41 

Note that Eq. (|5.4p has the same form as Ea. (|4.18|) therefore we can again treat it as a differential 
inequality as in the proof of Theorem 14.61 We obtain that 

A further application of Lemma 15.41 implies a similar result for the L^-norm of (p — po) a-nd the 
proof is finished. □ 

Remark 5.6. Here as in the case m = 2, when L is large, c' oc L^^ and we have the large L scaling 
of the rate proportional to in agreement with a perturbative analysis. However in 

this case, our arguments do not provide agreement with the constant of proportionality. We also 
note that by Theorem 12.41 we have that p{-^t) is uniformly Holder continuous in space and time 
for all t > T, where the Holder coefficient and exponent depends on 6, L and V. Thus we can 
bound, the L°°-norm of p — po by some power of its L^-norm. Hence the exponential convergence 
of Hp — Po\\l^{t^) implies the exponential convergence of ||p — PoIIl°o(t^)- However a bound along 
these lines is "even more" non-optimal since the two norms should, presumably, differ by a factor 
of L'^. 



6 Appendix: Proof of Lemma 12.51 

Proof of Lemma \2. 5\ We will do the comparison between p~ and w first; the comparison between 
p+ and w can be done in the same way. 

First note that w also satisfies Eq. p.Qp with $ = 0, therefore the inequality Eq. (|2.25p also hold 
for w, namely 

w- p- > -Cia, (6.1) 

where Ci depends on m, d. Mi, M2. 

We define f := w — p" , and our goal is to obtain an upper bound for /. More precisely, we want 
to show there exists some constant C and /3 depending on m, d, Mi, M2, such that /(x, t) < Ca^ in 

n X [1,2]. 

Our strategy is as following. First, we claim that 

g{T) := sup f f{x,T)dx < C^a for all T e [0,2], (6.2) 

yiEn-J B(y,i)nn 

where Co depends on m, d, Mi,A/2. We will prove this claim momentarily. Once we have the 
claim, we know the space integral of f{x,t) in any unit ball is of order a, for < t < 2. To get 
f{x,t) < 0{a^) for t E [1,2], it suffices to show / is Holder continuous in space with exponent and 
constant that are uniform in time for all t G [1, 2], which is indeed true, since Theorem 11.2 of |DGV] 
guarantees this uniform Holder continuity of p~ and w for t € [1,2]. 

Now it suffices to prove our claim. It is proved by writing both equations in weak form, choosing 
an appropriate test function and applying the Gronwall inequality. By writing both Eq. p.l9|) and 
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Ea. (|2.26p in weak form and subtracting the two equations, we arrive at 

/ f{x,T)(p{x)dx = [ f{x,0)ip{x)dx+ f f {w"" - p-\p-r^ -^)A^{x) + Ma^{x) dxdt, (6.3) 
Jh Jn Jo Jn 

^ ^ ^ ^ ^ ^ ^ ^ ^ 

where (p E C^{(l) is a test function chosen as follows. For a fixed T > 0, there exists z Q Cl, such 
that the maximum of /gj-y j^-jp^ f{x,T)dx is achieved at z. We then define 

(p{x) :^ fi * h^{x), 

where is a standard moUifier supported in i?(0, j^), and 

{l-|a;-z|V2 for|x-z|<l 
2)2/2 forl<|x-z|<2 (6.4) 
for |x - z| > 2 

For such ip, we have < (p < 1, inside the ball B{z, 1) and (pdx < \B{z, 3)| < 6''. 

To estimate Ji, note that (p{x) > 1/3 in B{z, 1), and f{x, T) + Cia > in 17, which implies 

7i = / [f[x,T) + Cia)Lp{x)dx — I Ciaip{x)dx 
Jn Jn 

> \ I {fix,T) + Cia)dx-6'^Cia 



3 



B{zA)nn 



3 

For I2, since /(x,0) = {:^)"^' p{x,Oy-'"' - p(x,0), we would obtain f{x,0) < Caa, where C2 
depends on m, ||/5(-, 0)||oo and c, (hence depends on m, d, Mi, M2), which yields 

I2 < Caa / ip{x)dx < 6'^C2a. 



Now we start to estimate I3. Due to the definition of in Eq. (j2.15p . we have m~ — m < 
2(m — l)ca. Also, we can derive some a priori bound of p~(x,t) and w{x,t) for t G [1,2], which 
depend on to, d, Mi, M2. Then we have 



w^^-p-\p-r'-' 



<C:i\w - p-\ + Cia inf^x[0,2], 



where (73,6*4 depends on m,d, Mi, M2. Together with the fact that \A(p\ is bounded, in particular 
by d, in B{z, 3) and vanishes outside of B(z, 3), we obtain the following bound for I3: 



h < {C3\f\ + Cia)\Ap\+ Maipdx 



n 



< dCs \f{x,t)\dx + 6''{dCi + M)a 

JB{z,3)nn 



1=1 JBiz+Xi,!) 



< dCsyi \f{x,t)\dx + 6'^{dC4 + M)a, (6.5) 

i)nn 
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where in the last inequahty wc denote by Cd the number such that B{Q, 3) can be covered by Cd 
numbers of balls of radius 1, centered at xi,...a;c^ respectively. Note that Cd is a constant only 
depending on d. 

Finally, we wish to control /g(^_,.^. ^j^q \ f\dx. Note that / > — Cia implies |/| < / + 2Cia, which 
yields 

/ \f{x,t)\dx < [ fdx + 2'^2Cxa 

J B{z+Xij)nn J B(z+xi.i)nn 

Plugging it into Ea. (|6.5p . wc obtain 

/3 < dCsCdgit) + {dCsCd2'^+^Ci + &'dC4 + Q'^M)a 
By putting estimates of /i, /2, ^3 together, we have 

9{T) <C5 [ g{t)dt + C&a for T e [0, 2] 
Jo 

where C5, Ce only depend on m, d, Mi, Af2- And for initial data, we have .g(0) < \B{0, 1)| sup^ f{x, 0) < 
2'^C2a. By Gronwall inequality, we have g{T) < Cqu for all T e [0,2], where Co only depends on 
m, d, All, M2, hence our claim Eq. (|6.2p is proved. □ 

Acknowledgements 

L.Chaycs was partially supported by NSF grant DMS-0805486, I.Kim was partially supported by 
NSF grant DMS-0970072, and Y.Yao was partially supported by NSF grant DMS-0805486 and 
DMS-0970072. 

References 

[BRB] J. Bedrossian, N. Rodriguez and A. Bertozzi, Local and global well-posedness for aggre- 
gation equations and Patlak-Keller-Segel models with degenerate diffusion, Nonlinearity, 
24(2011): 1683-1714. 

[Dib] E. DiBencdetto, Continuity of weak solutions to a general porous media equation, In- 

diana Math. Univ. J., 32(1983): 83-118. 

[DGV] E. DiBenedetto, U. Gianazza, V. Vespri, Harnack estimates for quasi-linear degenerate 
paraboUc differential equations. Acta Math., 200(2008), 181-209. 

[BS] A. Bertozzi and D. Slepcev, Existence and uniqueness of solutions to an aggregation 

equation with degenerate diffusion, Commun. Pure Appl. Anal, 9(2010): 1617-1637. 

[BH] M. Bertsch and D. Hilhorst, A density dependent diffusion equation in population dy- 

namics: stabilization to equilibrium, SIAM J. Math. Anal., 17(1986): 863-883. 

[BCM] Silvia Boi, Vincenzo Capasso and Daniela Morale, Modeling the aggregative behavior 
of ants of the species polyergus rufescens. Nonlinear Anal. Real World Appl, 1(2000): 
163-176. 



22 



[BCM2] Martin Burger, Vinccnzo Capasso, and Danicla Morale, On an aggregation model with 
long and short range interactions. Nonlinear Anal. Real World AppL, 8(2007): 939-958. 

[BF] M. Burger and M. D. Francesco, Large time behavior of nonlocal aggregation models 

with nonhnear diffusion, Netw. Heterog. Media, 3(2008): 749-785. 

[CJMTU] J. A. Carrillo, A. Jiingel, P.A. Markowich, G. Toscani, A. Unterreiter, Entropy dissipa- 
tion methods for degenerate parabolic problems and generalized Sobolev inequalities, 
Monatsh. Math., 133 (2001): 1-82. 

[CP] L. Chayes and V. Panferov, The McKean-Vlasov equation in finite volume, J. Statist. 

Phys., 138(2010): 351-380. 

[HV] M. A. Herrero and J. L. Velazquez, Chemotactic collapse for the KcUer-Segel model, J. 

Math. Biol, 35(1996):177-194. 

[KS] G. Karch and K. Suzuki, Blow-up versus global existence of solutions to aggregation 

equations, Appl. Math. (Warsaw), 38(2011): 243-258. 

[K] I. C. Kim, Erratum: "Degenerate diffusion with a drift potential: a viscosity solutions 

approach", Discrete Contin. Dyn. Syst., 30(2011): 375-377. 

[KL] I. G. Kim and H. K. Lei, Degenerate diffusion with a drift potential: A viscosity solutions 

approach, Discrete Contin. Dyn. Syst., 27(2010): 767-786. 

[KY] I. Kim and Y. Yao, The Patlak-Keller-Segel Model and Its Variations: Properties of 

Solutions via Maximum Principle, SIAM J. Math. Anal., 44(2012): 568-602. 

[LSU] O. A. Ladyzenskaja, V. A. Solonnikov, N. N. Uralceva. Linear and Quasilinear Equa- 

tions of Parabolic Type. AMS, Providence, RI (1968). 

[O] F Otto, The geometry of dissipative evolution equations: the porous medium equation, 

Comm. PDE., 26(2001): 101-174. 

[S] Y. Sugiyama, Global existence in sub-critical cases and finite time blow-up in super- 

critical cases to degenerate Keller-Segel systems, Diff. Int. Eqns., 19(2006):841-876. 

[TBL] Ghad M. Topaz, Andrea L. Bertozzi, and Mark A. Lewis, A nonlocal continuum model 

for biological aggregation. Bull. Math. Biol, 68(2006): 1601-1623. 

[V] J. Vazquez, The Porous Medium Equation: Mathematical Theory, Oxford University 

Press, 2007. 



23 



